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Elementary excitations in the multipole ordered state, which models the phase III in CeBe, are 
investigated by means of a generalized Holstein-Primakoff formalism. When different kinds of 
nearest-neighbor exchange interactions between multipoles are comparable to each other, orbital- 
flip excitations exhibit almost one-dimensional dispersion along the z axis. With high symmetry 
of the interactions, zero modes appear due to the macroscopic degeneracy of the ground state. 
The next-nearest-neighbor dipole-dipole interaction, which stabilizes the magnetic order of the 
phase III in CeBg, lifts the degeneracy and leads to gapfuU excitation spectrum. When the 
octupole-octupole next-nearest-neighbor interaction exists simultaneously, the spectrum shows 
softening at F and Z points. These excitations may be probed by neutron scattering and ultra- 
sonic measurements. 
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§1. Introduction 

Among many fascinating compounds with orbital de- 
generacy, CeBg shows an interesting entanglement of 
dipole, quadrupole as well as octupole moments. In 
particular, the entanglement gives a natural explanation 
for the longstanding inconsistency between the neutron 
diffraction and the B^^-NMR measurements of the phase 
11.^-' Namely, the field-induced octupole moment of T2- 
type, which arises through the mode mixing with the 
dipole moment, causes an extra internal field on the B 
nucleus sites. It explains the unexpected field-angle de- 
pendence in the NMR measurement. 

In discussing the transition from the antiferro- 
quadrupolar (AFQ) phase, called II, to the magneti- 
cally ordered phase III, the mode mixing between dipoles 
and octupoles also plays an important role.^' The 
present authors have made systematic analysis of multi- 
pole modes under the AFQ ordering. Then they have 
discussed effects of multipolar interactions on magnetic 
phases within the mean-field approximation. In their 
model Hamiltonian, the pseudo-dipole-type interactions 
for the next-nearest neighbors (NNN) are indispensable 
to account for the observed non-coUinear magnetic order- 
ing. The usual exchange-type interactions for the nearest 
neighbors (NN) are also important. 

What they found are the foUowings:^^-' (i) the com- 
petition between the NN dipole-dipole and octupole- 
octupole interactions strongly suppresses transition tem- 
peratures for simple AF magnetic order, (ii) the ferro- 
magnetic NNN dipole-dipole interaction plays a domi- 
nant role in stabilizing the non-coUincar ordering in the 
phase III, and (iii) the AF NNN octupole-octupole in- 
teraction is expected to be as strong as the dipole-dipole 
interaction since the phase III' is observed for relatively 
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weak magnetic field. 

For magnitude of the coupling constants in the model 
Hamiltonian, a reliable first-principle calculation is not 
available at present. However, the nature of the interac- 
tion parameters should be reflected in characteristics of 
elementary excitations. The purpose of the present pa- 
per is to derive spectra of all elementary excitations and 
to investigate how they reflect the details of interactions. 
The analysis also provides us with a typical example of 
collective excitations for systems where both the spin 
and the orbital degrees of freedom are involved. ^'^^ In 
this respect, a comparison with orbital ordered mangan- 
ites"'^'*^^^^ would be useful, in which weak coupling be- 
tween spin and orbital degrees of freedom gives rise to 
independent spin and orbital waves. In the present case, 
the excitations all involve magnetic moments due to the 
strong spin-orbit coupling of / electrons. Therefore all 
excitations can be measured by neutron scattering mea- 
surement. Some low-lying modes may also couple with 
ultrasonic waves. An experimental study for CeBg"'^*-* 
indicates strange structure of low-lying excitations. De- 
tailed experimental studies would be helpful to determine 
interaction parameters in the model Hamiltonian. 

This paper is organized as follows. In §2 we set up the 
model Hamiltonian for CeBg.^^ Then, a generalized 

Holstein-Primakoff formalism is developed by taking the 
reference state as the mean-field ground state. In §3 
the dispersion relations of spin-orbital wave excitations 
are given for some typical parameters in question. The 
summary is given in the final section. Details of the 
Holstein-Primakoff transformation are explained in the 
Appendix. 

§2. Model and Formalism 

2.1 Model Hamiltonian 

The four states in the Fg level, which is the crystalline- 
electric-field ground state of Ce'^'*' ion, are represented in 
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terms of the basis \ Jz) of J = 5/2 as follows: 



|1±) 



4 



|2±) = 



±- 



(2.1) 



The entanglement of spin and orbital degrees of freedoms 
can be described most physically by multipole operators 
in accordance with the point-group symmetry.^) The 
multipole moments are expressed in a concise way by 
means of two Pauli matrices, cr and r. Here r acts on 
the orbital partners, while a on the Kramers pairs, and 
is called spin hereafter. To specify a multipole operator, 
an irreducible representation F with possible multiplic- 
ity, and a component 7 are required. The set (F, 7) is 
abbreviated as A. In the irreducible representation F, 
the subscript u represents the odd property under the 
time reversal, and g the even one. The explicit form of 
the operators are summarized in Table I, where we 
have introduced linear combinations of and as 



(2.2) 

Experimentally, the Fsg-type quadrupolar ordering 
(phase II) first occurs at Q = (1/2, 1/2, 1/2) in units of 
27r/a, where a is the lattice constant. As temperature de- 
creases, the non-collinear magnetic ordering (phase III) 
is realized. The magnetic ordering is coexistent with 
the AFQ ordering, which is characterized by the dou- 
ble k modulation vectors, i.e., fei = (1/4,1/4,1/2) and 
fe2 = (l/4,-l/4,l/2).2.3) 

This complicated magnetic structure can be under- 
stood by the fact that a lowering of the point-group 
symmetry brought about by the AFQ ordering mixes 
otherwise independent multipole modes. For exam- 
ple, in the case of Oxy order the following mode mixings 
arise: 



(i) r2„-F4„i(z), 

(ii) ^4u2{x,y)-r5uix,y), 



(2.3) 
(2.4) 



where F4tii and r4ti2 represent F4„ with multiplicity in- 
dex 1 or 2. Since the order parameters in (ii) are con- 
cerned with realization of the phase III, it is natural to 
expect that the r4u2-type dipole-dipole and the Fs^-type 
octupole-octupole interactions are the most important. 

We assume that the intersite interactions between i 
and j are classified according to the representations of 
the cubic group. The conduction electrons that give rise 
to the Kondo effect are not treated explicitly. ^°"^^) Then 
the Fourier transform of the NN interaction for the 3D 
simple cubic lattice is given by 



~JrS^^'{cosqx +cosqy + cosqz)/3. (2.5) 



To realize the AFQ ordering of Fsg-type, J^g should be 
positive and the largest among all interactions, and hence 
we dcinote = Tq . The strength of the NN interactions 
mediated by conduction electrons was discussed''*' •^3) 
on the basis of the group-theoretical argument, where 
it was concluded that there are only two independent 
couplings. 

In order to account for the observed non-collinear spin 
orientation, we introduce the NNN interaction K of the 



pseudo-dipole type:^^'^^) 

Kf;^^' = -Kr{S^^' - 3n7,.n7;)/12, (2.6) 

where riij is the unit vector across i and j sites. 
The pseudo-dipole interaction is considered for three- 
dimensional odd representations, i.e., F4„i, F4„2 and 
r^u- Note that the positive (negative) Kr favors (anti) 
parallel alignment along riij. In the q space this inter- 
action is written as 




(2.7) 



where 



2 cos{qy) cos{qz) - cos{qz) cos{qx 



■ cos{qx) cos(gy) 



Jxy 



-sin(g^)sin(gy), (2.^ 



and other components are given by cyclic rotation of q^ , 
qy, and q^. The microscopic origin of the pseudo-dipole 
type interactions has been discussed by taking into ac- 
count the RKKY or the superexchange interactions to- 
gether with the d-f exchange or the Hund's-rule cou- 

pUllg.24-26) 

Consequently, the Hamiltonian used in this paper is 
given by 



i^j AB 



q AB 



(2.9) 



with the interaction 
Dr(9) = 



Jr(9), (F === 2^i, 3g, 5g), 

Jr(q) + Kr(q), (F = 4ul, 4u2, 5u). 

(2.10) 

According to the mean-field analysis,^ the transi- 
tion temperature of the simple AF ordering is suppressed 
by the competitions between J2„ and J4„i and between 
J4„2 and Jsij. In the case of the perfect competition, 
i.e., Jiu = •^4ui and J4„2 = •/su, the observed non- 
collinear ordering is stabilized for |/^5u|//^4u2 < 1 with 
^4«2 > 0.^^) The corresponding order parameters at 
T = are given by 

(Xf2)j^,jp = yl^^^(^^^_1^0)7A/2, 

(Xf)MF = A_Mi(6, 1,0)7^2, 



5gz\ 



/MF 



(2.11) 

(V3 ± l)/2 and (a,6,c)* 



where we have defined A±_ 
denotes a column vector. The site dependences are ex- 
pressed by 



ii = exp(iQ • Ri 



(2.12) 



1 TT TT 

^ii = —j= [cos(A;i ■ Ri + -)+ cos{k[ ■ Ri - -) 

- cos(A;2 ■Ri + '^)+ cos(fe2 • Ri - ^)] , (2.13) 
which specify the quadrupole sublattice {^i = ±1) and 
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the spin direction {fii = ±1) along the easy axis of each 
quadnipolc sublatticc. Here we have defined = Q — 
kn {n = 1,2). The real-space arrangement of {^i^fJ^i) is 
shown in Fig. 1. These mean- field solutions can be used 
as the reference state of elementary excitations in the 
Holstein-Primakoff formalism as will be developed in the 
next subsection. 

2.2 Generalized Holstein-Primakoff transformation 

For the purpose of investigating the elementary excita- 
tions, we first develop a generalized Holstein-Primakoff 
transformation for the Hamiltonian, (2.9). 

Let us denote the basis state for each site i by \n), 
where the quantization axes are common with all sites. 
Note that the multipole operators listed in Table I have 
been defined in this coordinate system. In the presence 
of a long-range order, especially a non-coUinear one, it is 
useful to introduce a local coordinate \n)i, which is taken 
according to be the ordering direction of the classical 
multipole. The local coordinate has been indicated by 
the subscript i of the ket vector. Practically, this can be 
done by the site-dependent unitary transformation, 



(2.14) 



In the case of the AF Heisenberg model for instance, 
such unitary transformation is given by the rotation tt 
with respect to the x axis on one of two sublattices. The 
unitary transformation is defined as such that the resul- 
tant local state \n)i becomes an eigenstate of the local 
mean-field Hamiltonian, 



MF 



(2.15) 



and its eigenvalues e„ are common to all sites. The 
averages of the physical operators with respect to the 
ground state reproduce the mean- field solutions, i.e., 
(0|X,^|0), = (X,^)mf. 

We introduce the Holstein-Primakoff bosons, b'^^, 
which describe local excitations from the ground state 
|0)i to excited states \n)i. Since we work with the local 
coordinate, the structure of the local excitations is com- 
mon for all sites. In terms of these bosons, the Hamilto- 
nian (2.9) is expressed (see Appendix for details) up to 
the second order as 



W = /f--iVM°°o"(0) 



/ , / .'''■MM 
m M 

■Mi-k) 6|„,M(fe)) X 



E' 

k 



(k) 



mn MN 



^Whik) 
^Mwik) 



^Mwik) 
^Mlkik) 



bl-Ni-k) 
bn,N{,k) 



(2.16) 



where the summation over k is restricted to the magnetic 
Brillouin zone, i.e., k is regarded as the reduced wave 
vector. Here we have defined the bosons in the k space 
with the band suffix M as 



bl.Mi'') 



m,k+GM ' 



where the reciprocal lattice vectors Gm are defined as 

such that k -\- Gm spans the original Brillouin zone, and 
G{) — 0. The matrix elements of ^l and A arc given by 

0^"^(fc) = AO:W(0) + MS™"(fc) - <5„„MoXm(0), 
= M~(fe), (2.18) 

with 

^M^Nik) = EE-^fe+Gif (^G„-Gjf)mn(^Gic-Gjv)fe'- 
AB K 

(2.19) 

Note that the ground state energy except for the zero- 
point energy has been subtracted from the Hamiltonian. 
By the following Bogoliubov transformation: 

as(fe) =EE["mM(fe)&m,-M(-fe) + <M(fe)CM(fe)]. 

(2.20) 



m M 

the Hamiltonian can be diagonalized as 



^ = E ' E ^^(^) Uc^Mk) + I) . (2.21 



The dispersion relation ijJs{k) and the coefficients and 
are obtained by the eigenvalue equation, 



\uj-F^lik)\=0, 



(2.22) 



^MivW - A^„^(fe) -^ZMk) ) ■ ^^■^'^> 

2.3 Local coordinate in the case of CcBq 

In order to rotate to the ordering direction, let us in- 
troduce the following unitary transformation, 

[/' = exp(-i7r/XiC7f/4)exp(-i7r^iC7f/8)(5exp(-i7r^iTf/4), 

(2.24) 



with 



Q 



( 1 





V 



\ 

1 


) 



(2.25) 



As will be shown shortly, exponential operators rotate 
spin and orbital states so that contains diagonal 

components and Tf only. The matrix Q interchanges 
the second and forth spin-orbital states. Then, the mul- 
tipole operators are transformed to 



In 



3g 



xT = 



i(x/3a+ - a^rn, |(^/3a^ - a+rf ), -^rf 



(2.26) 



(2.17) 

where we have defined af = {af ± iii^ia-)/\/2. 
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We factorize the mean fields in eq. (2.15) into its order 
parameter, eq. (2.11), and coupling strength, Ar as 

^4„2 ^ (Xf 2)^pA4„2, (2.27) 

</.f" = (Xf")MFA5„, (2.28) 

^5,z ^ (xf«^)MFA5,. (2.29) 

With = jj^ = 1, the mean- field Hamiltonian is diago- 
nalized as 

-"i — ^ a^XAu2-A_ ^ fTjA5„-TjA55 

(2.30) 

which depends neither on the sign of nor that of Hi 

as desired. The cigcnstatcs arc specified by |tz, crz)i, i.e., 
|0)i = !+,+)„ = \2)i = and |3)i = 

| — , — )i, and their energies are given by 

e+± = T{A\\iu2 + A^_\5u) - Ass, (2.31) 
e_± = tA+A_ (A4„2 - A5„) + Asg. (2.32) 

The schematic picture of local excitations is shown in 
Fig. 2, in which the oval represents the Oxy quadrupole 
and the thin (thick) arrow does the r4„2-type dipole (the 
Tsu-type octupole). Note that the dipole and the oc- 
tupole moments change their magnitude in the excita- 
tion processes with orbital flip (n = 2,3), while they do 
not in the spin-flip process (n = 1). 

By using the matrix elements of multipole operators, 
(2.26) for eqs. (2.19) and (2.18), the dispersion relation 
is calculated by eq. (2.22). 

§3. Spin-Orbital Wave Excitations 

In the phase III, the magnetic unit cell has the tetrag- 
onal symmetry and is 16 times larger than the para- 
magnetic one. Thus, there are 16 x 3 = 48 modes in 
the magnetic Brillouin zone. The definition of the high- 
symmetry points of the magnetic Brillouin zone is given 
in Fig. 3. 

Let us first investigate the case of the NN interactions 
only, i.e., K-p = 0. Figure 4(a) shows the case of the 
SU(4) limit, i.e., Jp = Tq for all F. The dispersion rela- 
tions are completely one dimensional; u) varies only along 
the z axis.^'^-' There arc 16 zero modes, which involve 
only the spin-flip (n = 1) excitations. Since without the 
orbital flip, the dipole-dipole and the octupole-octupole 
interactions cancel with each other, the ground state 
has a macroscopic degeneracy, 2^, which constitutes the 
zero modes. Namely, the assumed non-coUinear ordering 
is never stabilized in the case of the perfect competition. 
On the other hand, the excitations with n = 2, 3 bosons 
are always accompanied with the orbital flip, and the 
corresponding collective modes have the finite energy. 

When we go away from the SU(4) limit as shown in 
Fig. 4(b), some modes acquire dispersion in the F-X- 
M plane. However, overall structure of excitations does 
not change. There still remains almost one-dimensional 
excitations. When we set Jzg = Jau2 = Jsu = and leave 
only the competition between J2u and J^ui, the bunch 
of excitations around co/Tq ~ 1 and 2 in Fig. 4(b) merge 
roughly into one with co/Tq ~ 2 as shown in Fig. 4(c). 
If we leave the other competition, J4„2 = Jsu = 0.9Tq, 



two bunches of modes in Fig. 4(b) split into four bunches 
along Z-A-R-Z axis {oj/Tq 1,1.5,2 and 2.5 as shown 
in Fig. 4(d)). 

In any case, orbital-flip excitations with u) being larger 
than Tq show almost onc-dimcnsional dispersion rela- 
tions. The spin-flip excitations constitute all zero modes 
in consequence of the perfect competition between the 
NN interactions and of the lack of the NNN interactions. 
If we discuss possible lifting of the macroscopic degen- 
eracy without the NNN interactions, we should include 
quantum fluctuations in the ground state. Such consid- 
eration is out of the scope of this study. 

Next, we consider the effect of the NNN interactions 
with the NN interactions being fixed close to the SU(4) 
limit, i.e., J2u = 'hg = -hui = ■Jiu2 = Jbu = O.gTg. As 
Kiu2 is switched on and increased, the zero modes go 
up without losing dispersionless character and they form 
Ising-like excitations with a gap of the order of 2X4^2 ■ 
The case of Kiu2 = 0.5Tq is shown in Fig. 5(a). The 
Ising-like excitations can be understood in terms of the 
strong anisotropy with respect to the rotations of the 
dipole and the octupole moments in the presence of the 
quadrupolar ordering. 

In order to explain the stability of the phase III' 
in magnetic field, it has been proposed that the AF 
octupole-octupole interaction, K^^, should be as strong 
as the dipole-dipole interaction, i.e., K^^ ~ — -f^4u2-^^' 
Such situation in the coupling constants suggests that 
the Fsu-type octupole is a reasonable candidate for the 
order parameter in the phase IV, which was found re- 
cently in Cea;Lai_a;B6 with x ~ 0.75.^^"^^^ In the case 
of Kiu2 = 0.5Tq and K^u = —QA^Tq, the Ising-like ex- 
citations become quite dispersive and show softening at 
F and Z points as shown in Fig. 5(b). If one increases 
|-^5ii| beyond i^4„2, the energy of the soft modes becomes 
complex. This means that the assumed ground state is 
unstable against exponentially growing modes. 

§4. Summary 

We have investigated the spin-orbital wave excitations 
for the phase III of CeBg. Using the mean- filed solution 
as the classical ground state, we calculate the dispersion 
relations in a generalized Holstein-Primakoff formalism. 

When the nearest-neighbor dipole-dipole and octupole- 
octupole interactions are comparable, the dispersion of 
orbital-flip excitations is almost one-dimensional along 
the z axis. The spin-flip excitations constitute zero 
modes owing to the macroscopic degeneracy of the 
ground state. 

As the next-nearest-neighbor dipole-dipole interaction 
which stabilizes the magnetic order of the phase III in- 
creases, the; zero modes go up keeping the dispersion- 
less character and form gapfuU Ising-like excitation spec- 
trum. When the dipole-dipole and the octupole-octupole 
interactions exist simultaneously, the Ising-like excita- 
tion shows softening at F and Z points. We expect 
that these modes can be probed experimentally. De- 
tailed measurements of excitations would reveal the na- 
ture of spin-orbital excitations and relevant intersite in- 
teractions. 
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Appendix: Details of Generalized Holstein- 
PrimakofF Transformation 

In this appendix, detailed derivation of (2.16) is given 
by introducing a generalized Holstein-Primakoff trans- 
formation for the Hamiltonian, (2.9). 

Prom now on we work with the local coordinate which 
is defined by (2.14). By introducing the Hubbard oper- 
ator at i site, 51"" = \m){n\i, the Hamiltonian, (2.9), is 
expressed as 

H=^J2Y1 "Sf , (A-1) 

ij mnkl 

where the intersite interactions are given in terms of the 
matrix elements of the multipole operators. In the local 
coordinate we obtain 



of MJ^"*^' is given by 



M™" 



kl 



D 



AB 



iXf)rr.Mf)^. (A-2) 



The generalized Holstein-Primakoff transformation is 
defined as (for m, n ^ 0) 



^00 



V ' 



~ b\nibnii 



(A-3) 
(A-4) 

(A-5) 
(A-6) 



where the prime means that the summation excludes £ = 
0. The commutation relation for the Hubbard operators, 



(A-7) 



is satisfied if one uses boson commutation relations for 
bii and 6].. The boson creation operator bj- represents 
the local excitation from the ground state to the state 
i. The symbol P is the expectation value of S^° with 
respect to the vacuum of bosons. Then we can perform 
a systematic expansion of the Hamiltonian in power of 
1/P. After the expansion we set P = 1 for the physical 
system. 

Up to 0{1/P), the Hamiltonian is expressed as 



2 



iVM°7 + VNj:j:'iMro"blk + h.c.) 



EE' 

kp mn 



Qmn L t /, 



^{^k"pbm,-kbn,p +i^-C.) 



(A-8) 



where N is the number of sites and the Fourier transform 



Mk,p' =J2^^''(^k-cl)mn{X^.p)kl. (A-9) 



Here we have defined 



k,p 



—SmnM'. 



0,p-k + ^0,p-k 



k,p 



: M 



k,p 



(A-IO) 
(All) 



If D^^ is real, then we have the following relations 

In eq. (A-8), the first term denotes the energy of the 

classical ground state. The second term linear in bo- 
son operators must vanish, since otherwise the assumed 
classical ground state is no longer stabilized. 

To diagonalize the Hamiltonian (A-8) with respect to 
fc, we introduce the reciprocal lattice vectors Gm such 
that the wave vector k within the magnetic Brillouin 
zone plus Gm spans the original Brillouin zone. By not- 
ing that and A™^ are proportional to Sp^k+Gu '^itih 
M corresponding to the ordering vectors, we obtain the 
Hamiltonian of eq. (2.16). 
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Table I. The multipole operators in the Fg subspace. 
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Fig. 2. The schematic picture of local excitations. The ovals rep- 
resent the Oxy quadrupole, thin arrows the F4ii2-type dipole, 
and thick arrows the Fs^-type octupole moments. 
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Fig. 1. The real-space arrangement of the phase III in the (001) 
plane with the coordinate z. The symbol {^i,fii) characterize 
the quadrupole sublattice (^i) and the spin direction (fii) along 
the easy axis of each quadrupole sublattice. 



Fig. 3. The high-symmetry points in the magnetic Brillouin zone 
with tetragonal symmetry. 
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(a) Jr = Tq for all F (b) J2a = J3g = J4ul = Ma2 = ha = 0-9Tq 
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(c) = J4ui = 0.9TQ (d) J4u2 = J5u = 0.9Tq 
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Fig. 4. The collective excitations without the next-nearest- 
neighbor interactions. The mean-field solution corresponding 
to the observed non-coUinear ordering is used for the classical 
ground state. 
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(a) J2u = J3g = J4ul = J4u2 = ha = 0-9Tq 
K4u2 = 6.5Tq,K5u = 
4.0 
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(b) J2u = J3g = J4ul = Ma2 = ha = 0-9Tq, 
K4u2 = 0.5Tq,K5u = -0.45Tq 
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Fig. 5. The collective excitations for (a) K4u2 = 0.5Tq and (b) 
K4u2 = 0.5Tq, i^sti = — 0.45rQ with the nearest-neighbor in- 
teractions close to the SU{4) limit, J2u = J3g = Jiul = Jiu2 = 
J5„ = 0.9Tq. 



